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ABSTRACT 

We examine the basic physics of type la supernova (SNe la) light curves 
with a view toward interpreting the relations between peak luminosity, peak 
width, and late-time slope in terms of the properties of the underlying explosion 
models. We develop an analytic solution of the radiation transport problem in 
the dynamic diffusion regime and show that it reproduces bolometric light curves 
produced by more detailed calculations under the assumption of a constant 
extinction coefficient. This model is used to derive the thermal conditions in 
the interior of SNe la, to address the issue of time dependence in the modeling 
of SNe la light curves and spectra and to show that these are intrinsically 
time-dependent phenomena. The analytic model is then used to study the 
sensitivity of light curves to various properties of supernova explosions. We 
show that the dominant opacity arises from line transitions and discuss the 
nature of line opacities in expanding media and the appropriateness of various 
mean opacities used in light curve calculations. 



Subject headings: sup ernovae: general, cosmology:distance scale, radiative 
transfer 



^General Studies Group, Lawrence Livermore National Laboratory, Livermore CA 94550 USA 



- 2 - 



1. Introduction 

The bolometric light curve is the simplest and most direct manifestation of type la 
supernovae (hereinafter SNe la). For many years it had been assumed that all Type la 
supernovae (hereinafter SNe la) were identical explosions, with identical light curve shapes 
and peak luminosities (c.f. [Woosley fc Weaver 19861) . While evidence for this uniformity in 
the data was never terribly convincing, the use of SNe la as the primary "standard candles" 
for cosmological distance measurement provided a powerful incentive for assuming this 
homogeneity. This in turn lead naturally to a search for an explosion model which might 
produce identical displays from the diversity of progenitors supplied by stellar evolution. 

It became clear from the light curve's rapid evolution that a relatively low mass object 
was involved, one with a short radiative diffusion time ( [Arnett 1982|) . The result of this 
search is what one might call the present "standard model" , the thermonuclear incineration 
of a carbon-oxygen white dwarf at the Chandrasekhar mass (see [Woosley fc Weaver 1986 



Arnett 1996| for details of this search and a review of various models). The Chandrasekhar 



mass provides a point of convergent evolution for various progenitor systems, offering a 
natural explanation of the assumed uniformity of display. The high densities attained at 
the centers of these objects provides as well a mechanism (as ill-defined as it may be at 
present) for their ignition. 

With the coming of age of various supernova searches, ( [Hamuy et al. 1993| , |Pollas| 
1994 , Evans, Van Den Bergh, fc McClure 1989 , [Barbon et al. 1993|) there has recently 



been an explosion in the availability of high-quality data. It is now generally recognized 
that SNe la exhibit a variety of light curve shapes, peak luminosities, and maximum-light 
spectra. Perhaps most significant has been the discovery of various regularities in the light 
curve data, the most famous of which we will call the "Phillips Relation" (hereinafter PR) 
( [Phillips 19931 ); the brightest supernovae have the broadest light curve peaks. There is also 



significant evidence ( [Vacca fc Leibundgut 1996a ) that the decline at late times is more 



gradual in the brighter supernova. 

This additional information provides new clues to the nature of these explosions. We 
develop here, from first principles, a theoretical framework for examining the formation of 
SNa la light curves and extracting underlying properties of the explosions. While there 
have been a number of studies of the light curve problem in SNe la ( [Hofiich, Khokhlov 



fe Wheeler 1995[ and references therein, [Weaver, Axelrod, fc Woosley 19SU[ ), none has 



included a detailed examination of the physics of the opacity or radiation transport in these 
explosions. 



We divide the present investigations into four parts. In the first we derive an analytic 
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solution for the co-moving frame transfer equation in homologous expansion. We use this 
solution to examine the relative importance of various terms in the equations and the 
suitability of various approximations which have appeared in the literature. The analytic 
solution also provides an important check on the accuracy of our subsequent numerical 
solutions. 

In the second part we examine the nature of the opacity. SNe la differ significantly 
from other astrophysical objects in their composition; they are entirely constituted of heavy 
elements. This allows a number of interesting effects, present at a low level in all objects, 
to take on a dominant role. We find that the extreme complexity of the atomic physics 
allows for a transport of energy downward in frequency which has a profound effect upon 
the radiative diffusion time. This transport also allows for a process, akin to thermalization 
but not mediated by collisions, which can lead to a spectrum which appears thermal but 
which need have little to do with the gas temperature anywhere in the object. 

In the third part we present several light curves calculated by solving the time- 
dependent, multi-frequency transport equations. We use these detailed solutions to examine 
the behavior of the mean opacity, and find that the most natural a priori choice, the 
Rosseland mean, significantly underestimates the true flux mean by up to a factor of five. 
The flux mean appears to be, at least in LTE, remarkably constant both in time and in 
radius. 

Finally, we exploit the constancy of the flux mean to examine the effect of varying the 
bulk parameters of the explosion on the shape of the light curve. For present purposes, 
SNe la can be characterized by specifying density and composition as a function of ejecta 
velocity. We will ignore, at least to begin with, how this structure was brought about. We 
will also ignore any departures from spherical symmetry — as will be seen, the problem is 
vexing enough in one spatial dimension. We note however that Wang et al. ( |199(j| ) found 
no evidence for polarization in their study of three SNe la, nor has there been any other 
evidence from optical polarimetry for gross asphericity in SNe la explosions. 

The post-explosion dynamics, which determine the density and velocity structure, are 
quite simple. As in all strong point explosions, the expansion becomes homologous after 
a time which is short compared with the bolometric rise-time, with a velocity gradient 
everywhere equal to the reciprocal of the elapsed time. The density structure is quite 
smooth, with a density profile not very different from p oc exp(— 4f/fo) (with vq some 
typical velocity oc [E / My/"^). There is usually a small amount of high- velocity material at 
the surface in which the density drops more rapidly; as this material is largely transparent 
long before the supernova is observed we may ignore this detail. To a surprising degree, all 
explosion models to date have this relatively simple structure. Thus, the dynamics can be 
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specified by the explosion energy and the mass of the ejecta (in most models the progenitor 
is completely disrupted). 

The other defining attributes of the explosion are in the composition of the ejecta. 
Since it is radioactive nickel which leads to any optical display at all, the amount of ^^Ni 
and the depth to which it is buried in the ejecta will obviously affect the light curve. 
The composition also affects the opacity, obviously a determining parameter in a problem 
concerning the escape of radiation. 

We will thus define a supernova by its total mass, explosion energy, ^^Ni mass and 
opacity, possibly allowing for variations in the spatial distribution of the latter two (we will 
find that this is relatively unimportant). With a simple means of producing a light curve 
from these four parameters, one may turn the problem around and use the light curve 
model to extract values for these parameters from observations of SNe la. 

We find that, with the exception of the total mass, variations in any of these basic 
parameters lead to a behavior of the light curve which is in the opposite sense of the PR. 
Varying the total mass, however, can lead to a sequence of light curves in which the PR 
behavior is reproduced. This is by no means a proof that the full richness of SN la light 
curve behavior cannot be obtained from Chandrasekhar mass explosions. It does however 
suggest that the total mass of the explosion may be a natural and simple explanation for 
observations. 

This is the first paper of a series. The next paper ( [Pinto fc Eastman 1996| ) explores 
the lightcurve behavior of simple models for SNe la. Subsequent papers will systematically 
explore the light curve and spectrum properties of specific models for Type la supernovae. 

2. A Schematic Type la Supernova 

In this section we develop a simple analytic model for the thermal evolution and light 
curve of a SN la. This will prove useful for estimating physical conditions in the ejecta at 
various times after explosion, and for illustrating the effect which changes in opacity, mass, 
energy deposition, and explosion energy have on the bolometric light curve. 

The ejecta of SNe la form an opaque, expanding sphere into which energy is deposited 
by radioactive decay at an exponentially-declining rate. Because the sphere is initially 
so opaque, this energy is converted into kinetic energy of expansion on a hydrodynamic 
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timescaleQ At the earliest stages the ejecta is so optically thick that the time it takes 
radiation to diffuse out is much longer than the elapsed time. The luminosity is therefore 
initially quite small. As time passes, the ejecta become more dilute and the diffusion time 
drops below the (ever-increasing) elapsed time. Since the rate of energy input declines 
exponentially with time, there is a peak in the light curve as soon as the injected energy 
has an appreciable chance to escape conversion to kinetic energy — when the diffusion time 
becomes comparable to the elapsed time. While the fraction of deposited energy which 
escapes conversion will continue to increase, this is more than offset by the decreasing 
energy deposition rate. 

Shortly after this peak, there will be a considerable amount of radiation still trapped 
and diffusing outward. Since the energy deposition rate is so rapidly declining, the 
luminosity will, for a time, exceed the rate of deposition until the supernova empties itself of 
this excess stored energy. Finally, as the rate of energy deposition, now from cobalt decay, 
declines more slowly and the diffusion time becomes small, the luminosity becomes equal 
to the instantaneous deposition rate. There are thus two milestones in the light curve. The 
first occurs near peak when the luminosity first rises above the rate of energy deposition. 
The second occurs when the excess, stored energy is exhausted and the luminosity falls to 
equal the instantaneous deposition. The elapsed time and the rate of deposition are easily 
determined. The first is obvious and the second comes from the decay of ^^Ni to ^^Fe and 
the transport of 7-rays - fairly simple physics. Determining the diffusion time is a far more 
complex matter, and most of the difficulty in producing synthetic light curves and spectra 
arises from correctly characterizing the transport and escape of thermalized radiation. 



Arnett showed in two elegant papers ( |1980| ; |1982|) that the ideas expressed above could 



be demonstrated by a simple analytic model which accounts for the deposition and escape 
of radiation from the expanding ejecta. This model predicted a bolometric light curve 
which was generically in good agreement with observed Sn la behavior. Starting from the 
thermodynamics of the trapped radiation, he showed that the luminosity at peak was equal 
to the instantaneous energy deposition rate under the assumption of constant opacity, and 
thus the first milestone occurs near peak bolometric luminosity. A number of assumptions 
were made which, for a first attempt, were quite reasonable, but which rendered suspect the 
precise predictions for any particular model explosion. These included a constant density 
structure, a constant opacity in both space and time, and the requirement that the radial 
distribution of the energy deposition was identical to that of the thermal energy. Thus, 
while he could vary the expansion velocity and the total mass, the effects of varying the 



^For a point explosion like a SN la, the hydrodynamic timescale is comparable to the elapsed time. 
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structure of the ejecta and of a realistic energy deposition profile were beyond examination. 
In mathematical terms, Arnett's solution was an eigenfunction expansion from which only 
the fundamental mode is retained. We shall have more to say on this later. 



We take a no-frills approach similar to Arnett's (|1982|) , while relaxing some of his 
more limiting assumptions so as to be able to address additional questions such as how the 
density structure and distribution of radioactive isotopes are manifested in the bolometric 
light curve. 

We start by writing down the first two frequency-integrated moments of the time 
dependent, co-moving frame radiative transport equation in spherical geometry. The first 
moment equation, for the radiation energy density, can be written, correct to all terms 
0(v/c), as (cf. [Mihalas fc Mihalas 198^) 



Dt 

The second frequency- integrated moment, for the radiation momentum, is 



00 



IDF dP 3P-E 2 dv v \ ^ 1 , 

+ + + ^hr + -K= — / XuF^^du. (2) 



Dt dr r \dr r J c Jo 

Here E, F, and P are the zero, first, and second frequency-integrated moments of the 
radiation field: the energy density, the flux, and the (isotropic) radiation pressure. Xu is 
the extinction coefficient, and rji, is the volume emissivity. These are formidable equations 
to solve directly (cf. [Eastman fc Pinto 1993| among others) and our goal here is to obtain a 



simple and approximate solution. 

The first and most important approximation we will employ is that the expansion 
is homologous. As already described, SNe la are strong point explosions; homologous 
expansion will be an excellent approximation if the energy released by ^^Ni decay does not 
strongly affect the dynamics of the expansion. The energy available from ^^Ni decay is 
3 X lO^^erg g~^. This corresponds to the kinetic energy of a gram of material traveling at 
nearly 2500 km s~^, or, equivalently, a velocity increment of the same magnitude over the 
velocity initially imparted by the explosion. The significantly greater decay energy available 
from decay all the way to ^^Fe is less relevant as most of the ^^Co decay energy is emitted at 
times when the supernova is becoming optically thin. Since the observed expansion velocity 
of SNe la is in excess of 10^ km s~^, we expect that this additional source of energy will 
have a modest, but perhaps not completely negligible, effect upon the velocity structure. 
Furthermore, as the time to maximum light, tmax, is observed to be at least twice as large 
the ^^Ni decay time, most of the hydrodynamic effect of ^^Ni decay will have occured 
prior to a supernova's discovery. If we take the ejecta's density structure from an explosion 
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calculation which has allowed this additional energy to accelerate the ejecta for the first 
few days we will have taken this effect sufficiently into account. We will therefore take the 
outer edge of the supernova, or at least of a fiducial mass shell which contains virtually all 
of the mass, to be at a velocity Vmax and a radius 

R{t) = -Ro + Vmaxt (3) 

where Rq is the initial radius of the progenitor. For this type of expansion law there is an 
associated time scale 

tsc Rq/ ^max ■ (4) 

which will be one of the parameters of the solution. 

A major simplification is the so-called Eddington Approximation, wherein the radiation 
field is everywhere isotropic: E = 3P. This is certainly valid during the early, optically 
thick stages of evolution, but breaks down when the ejecta become transparent. The error 
that this assumption introduces at late times lies in the energy distribution, and has much 
less effect on the bolometric luminosity which is the only link this simple analytic model 
has to SN la observations. We note in passing that this assumption has dire consequences 
for the calculation of the energy deposition, where the deposition rate is proportional to the 
gamma-ray energy density. 

At the temperatures and densities of maximum light SNe la, the gas energy density 
is less than the radiation field energy density by a large factor. The radiation field energy 
density is 

/ T \^ 

perad - aT^ ~ 1210 — (5) 

which greatly exceeds both the thermal kinetic energy density 

. + ^ 0.4 [^^) (^) (I) (6, 

and the ionization energy density 

where Na is Avogadro's number, A is the mean mass per nucleon, i is the average ionization, 
and {E) is the mean ionization energy. 

The dominance of radiation over internal energy permits us to ignore the gas internal 
energy and set 

poo 

/ [Anrj^ - cXuE^] dv = e (8) 

JO 
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where e is the volume rate of 7-ray deposition. Equation | is equivalent to saying that 
as soon as high-energy radiation from decay is absorbed, it is immediately re-radiated as 
thermal emission. The mechanism by which this happens is collisional; 7-rays Compton 
scatter, producing high-energy electrons which then rapidly transfer their energy to the 
plasma. This occurs on a timescale which is short compared with any other timescale 
important to the problem of energy transport. 

Much of the radiation transport in the peak-light phase of the light curve occurs as 
diffusion. As discussed earlier, the diffusion time for radiation to escape the supernova 
is, at first, long compared with the hydrodynamical (elapsed) time. This puts us in the 
dynamic diffusion regime of ( [IVlihalas fc IVlihalas 1984| ) . Following the discussions referenced 
therein, it is important when in this regime to retain all of the terms in the radiation 
energy equation to 0(v/c) . However, in the frequency integrated radiation momentum 
equation, equation |], it is appropriate to discard all time- and velocity-dependent terms. 
This difference in treatment is intuitively evident when one realizes that we are vitally 
interested in determining the energy density and its flow within the supernova, yet the 
radiation momentum has little effect upon the supernova's dynamics after the first few 
days. The radiation momentum equation is thus reduced to the familiar diffusion form 

p__J^dE_ 
3x Or ' 



where x is an appropriately defined frequency-averaged mean opacity (see section 
below). Using this and the previous approximations in the radiation energy equation, the 
transport equation becomes 

In general, equation |10| is too complicated to solve analytically. However, for certain 
conditions specified below the spatial and temporal parts of the solution may be separated. 



and equation equation 10 reduced to two ordinary differential equations. 



It is convenient to homologously scale all of the remaining quantities in terms of the 
time t and x = r/R, the (dimensionless) fractional radius. Since the gas is radiation 
dominated, E oc R'"^ oc (1 + t/tsc)~^, we write 



E{x,t) = S{x,t)Eo 



-V 



Mt), 



EoHx)4>{t)(^] , (11) 



- 9 - 



where ip{x) describes the radial variation, the temporal variation, and Eq is the initial 
energy density. 

The density can be written as 

p(r,t)=poP(x)(^^y, (12) 

where p{x) is the radial profile of the density, normalized so that p(0) = 1. 

The extinction coefficient x is the mass opacity coefficient k times the density. We allow 
K to have an intrinsic radial dependence, described by Kj(^X^ J clS well as a time dependence. 

City. 

X{r,t) = KoK{x)C{t)pop{x) (^-^^ ■ (13) 

Separation of variables is possible only if the opacity does not depend upon the energy 
density (i.e. the temperature). Remarkably, conditions may actually conspire to produce 
a mean opacity which is roughly constant with both time and depth through the ejecta 
(section 



The volume energy deposition rate e will scale as 

.(,,.?^,,)A(,,(j|y (14) 

where eo = Em + Eco = 1-73 MeV + 3.69 MeV = 5.42 MeV is the total energy available 
from decay, per atom of ^^Ni, and A{x,t) is the dimensionless energy deposition function 
which results from 7-ray transport and escape. The total production rate of decay energy 
function of t is described by 6{t), given as 

e{t) = eo-i (Eme-'/^^' + (e"*/-^" - e'*/^-) I • (15) 

It is convenient to define the total energy available from the ^^Ni — i>^^Fe decay chain 
(excluding neutrinos) in terms of the total kinetic energy of the gas, 

3M„ieo 

and the initial diffusion time from the center as 

= , (17j 

c 

where xo = Po^^o- 
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Substituting equations |TT], |T2], 0, [T^ and |T5] into equation |T0| gives 



R{t) 



x"^ dx \x dx 



where we have combined the spatial shape of the density and opacities into the function 
x{x) = pk. Equation |18| is the principle equation which describes the evolution of the 
radiation field energy density. 

Next we must specify the boundary conditions. At the center there is a reflection 
boundary condition where the flux vanishes, or equivalently, ip'{0) = 0. At the surface we 
use a solution to the plane-parallel grey atmosphere problem, assuming that the thickness 
of the surface layers is small compared with the radius. This can be expressed as 



ip{x) 



At the outer edge, r = 0, and we have 



and thus 



The boundary condition is then 



ipix) = -ip{l) r + - 



2 dtp 



x=l 



dr 
dx 



3 dx 

Since the optical depth to the surface from radius r is 

pudr 



(19) 

(20) 
(21) 

(22) 



PoKoR{t) / p{x')K{x')dx' 



we have 



dx 



-poKoR{t)p{l)~K{l) 



(23) 



(24) 



x=l 



It is more convenient as a boundary condition to require the solution to go to zero at some 
radius beyond x = 1. Extrapolating equation ^linearly, we find that this is equivalent to 
requiring 'ip^Xo) = at 

xo = 1 - 2/3 



dr 



(25) 



1=1 
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The value of Xq increases as density declines and strictly speaking, this will introduce a 
time dependence into the spatial solution which violates the conditions making equation [18 
separable. However if we consider this to be a slow, quasi-static change, the solution 
obtained by separation of variables is not too inaccurate, and will be adequate for our 
needs. We will touch on this again. 



To solve equation 18 we follow the usual procedure for separation of variables and first 



find a solution to the homogenous equation, where energy deposition from decay is set to 



zero. In the absence of any sources, equation 18 can be written as 



'^^r,((d^^l(C'4)- (26) 



R{t) (J) x'^ip dx \ x dx ^ 

Since the terms on the left hand side of this equation depend only on t, while those on the 
right hand side depend only on x, each must be equal to a constant independent of either x 
or t. Let this constant be a. We can then write, for the spatial part. 



d ( x"^ dip' 



x'^ip{x) dx \x dx 
For X = 1) the solution can be written 



-a. (27) 



where the eigenvalue a depends upon the total optical depth. For large opacity, dr/dx is 



large, and the boundary condition equation 24 approaches the radiative-zero condition, 
?/'(l) = 0. For constant Xi the radiative- zero eigenvalues are a„ = ra^vr^ and the 
eigenfunctions are 

ijn{x) = V2^ ^. (29) 

X 

Given a„, the temporal part of the solution, (j)n{t) is determined by the homogenous 
equation 

^r,C(t)^ = (30) 

R{t) (pn 

For ({t) = 1 (opacity does not change with time), the solution can be written as 

,„,„.exp(-^^^(l^^)). (31) 

When X is an arbitrary function, there are no analytic solutions to equation 
and we must determine the eigenfunctions numerically. Eigenvalues are determined by a 
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Rayleigh-Ritz procedure, and a discrete representation of their corresponding eigenfunctions 
is obtained by relaxation. The basin of convergence to a desired eigenfunction for this 
process is surprisingly small; for most interesting p{x), eigenvalues must be determined to 
better than a percent for the resulting relaxation to converge to the desired solution. We 
prefer to normalize the solutions such that the functions ipn are orthonormal with respect 
to the inner product 

<f\g>= [' f{x)g{x)x^dx. (32) 







As the solution progresses in time, the spatial solution changes because the value of xq 



increases with decreasing optical depth. To avoid the need of a new solution of equation |27 
at each time, in much of the discussion which follows, we will take the radiative zero 
solution. This allows a single set of eigenvalues to be used at all times. For more realistic 
calculations, we note that the change in eigenvalues due to changes in x over time is slow 
and continuous. We may thus continuously re-solve the eigenvalue problem as we evolve 
the solution in time. 

As an example. Table |1| lists the first 25 eigenvalues for the radiative zero solution with 
X = e~^^, with = 1 to 4. A selection of these functions is shown in Figure 



Turning now to our original, inhomogeneous transport equation, equation |l^, the 
general solution for E{x,t) is an expansion in terms of the eigenfunctions ipn'- 

oo 

E{x,t) = J2 Em^m{x)(f)m{t)- (33) 
m=l 

If we substitute this expansion into the inhomogeneous equation |TB|, multiply by x'^ipnix) 
and integrate from x = to 1, we get 

0„ + — (1 + t/tsc) m-'(pn = m) (1 + t/tsc) < A|Vn > • (34) 

where < A|?/'„ > is the overlap integral of A{x,t) with eigenfunction ipn{x)- 

For simple forms of ({t) and < A\ipn >, equation ^ can be integrated analytically. It 
is, however, a straightforward matter to integrate this equation numerically, and one need 
not be limited by the analytic integrability of these two functions. 

We must now put back the dimensional constants to be able to calculate a light curve. 
Starting with the definition of the flux equation ^ we have for the luminosity 
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Table 1 
Eigenvalues for p{x) = e 



mode 


k=0 


k=0.25 


k=0.5 


k=l 


k=2 


k=4 


1 


9.8696o 


1.1878i 


I.425O1 


2.0298i 


3.9256i 


1.16372 


2 


3.9478i 


4.5567i 


5.2460i 


6.9007i 


1.15962 


2.91752 


3 


8.8827i 


1.01522 


1.15732 


1.49222 


2.40892 


5.64532 


4 


I.579I2 


1.79782 


2.04142 


2.61152 


4.14772 


9.39372 


5 


2.46742 


2.80372 


3.17732 


4.04882 


6.378I2 


I.4I8I3 


6 


3.553I2 


4.03282 


4.56512 


5.80432 


9.10122 


2.00153 


7 


4.83622 


5.48522 


6.205O2 


7.87822 


1.23173 


2.69OO3 


8 


6.31672 


7.I6IO2 


8.09692 


I.O27I3 


1.60273 


3.48373 


9 


7.99462 


9.06022 


I.O24I3 


1.29823 


2.O23O3 


4.38283 


10 


9.87OO2 


1.11833 


1.26373 


I.6OII3 


2.49273 


5.38743 


11 


1.19433 


1.35293 


1.52853 


I.936O3 


3.01173 


6.49743 


12 


1.42133 


I.6O983 


I.8I863 


2.30273 


3.58023 


7.713O3 


13 


I.668I3 


I.889I3 


2.13383 


2.70133 


4.I98O3 


9.03423 


14 


1.93463 


2.19073 


2.47433 


3.13173 


4.86533 


1.04614 


15 


2.22O83 


2.51473 


2.84OO3 


3.594I3 


5.58193 


1.19934 


16 


2.52693 


2.86IO3 


3.23093 


4.08833 


6.348O3 


1.36324 


17 


2.85263 


3.22973 


3.647I3 


4.6I443 


7.16353 


1.53754 


18 


3.I98I3 


3.62073 


4.08853 


5.17243 


8.02843 


1.72254 


19 


3.56343 


4.O34I3 


4.555I3 


5.76243 


8.94283 


I.9I8O4 


20 


3.94843 


4.46993 


5.04693 


6.38423 


9.9O663 


2.I24I4 


21 


4.35323 


4.92793 


5.56403 


7.03793 


I.O92O4 


2.34O84 


22 


4.77783 


5.40843 


6.10633 


7.72353 


1.19834 


2.568O4 


23 


5.222I3 


5.91123 


6.67393 


8.44IO3 


1.30954 


2.80594 


24 


5.68623 


6.43643 


7.26673 


9.19053 


1.42564 


3.05434 


25 


6.I7OO3 


6.98403 


7.88483 


9.97193 


1.54674 


3.31334 


ratio 


1.0 


1.1322 


1.2782 


1.6163 


2.5063 


5.3643 



Table 1: The first 25 eigenvalues for the spatial equation ^ The final row is the ratio of subsequent 
eigenvalues to n^7r^, the asymptotic limit. 
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Fig. 1. — Eigenfunctions for mode numbers 1, 2, 5, and 10 for exponential density laws with constant 
if = 0,4 
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Because the boundary conditions are ipi^) = and ip'^Xo) = 0, there is no scale to the 
problem, and we are free to impose a third condition on the overall normalization of the 
solution. We can compare the luminosity with the energy deposited at late times by noting 
that when the timescale over which the energy deposition 6{t)A{t, x) changes becomes long 
compared with the diffusion time, the go asymptotically to 

0„ = ^e^(t)C(t)(A|V'n). (36) 
an 

If we integrate equation over volume, we find that 



dip \ 

= -a„4(l), (37) 



x(x) dx ^ 



x=l 



where the total radiation energy interior to x is 



X 



J„(x) = X ijjndx. (38) 







Putting these two results into the expression for the luminosity and using the definitions of 
Td and e gives 

oo 

L{t) = 3M^,eo9{t) J2 0n(t)(A|V^n)/n(l) (39) 



n=0 

If A{x,t) is constant in x, we can let A{x,t) ~ A(t), and the sum becomes 

1 2 



A(t)E 



X^lpndx 



n=0 







(40) 



Requiring, then, that the solutions ipn to be normalized such that this is A(i(:)/3 leads 
finally to 

L(t) =MW(t)A(t), (41) 
the instantaneous energy deposition. 

In order to examine the effect of including an increasing number of modes on the shape 
of the light curve it is convenient to renormalize the energy deposition such that the correct 
total amount of energy is deposited per unit time into whatever modes are included in the 
calculation. We therefore divide the energy deposition factor {A\ipn) by the quantity 



/ = (42) 

/q x'^Adx 



This has an aliasing effect of overestimating the power in the included modes just enough 
to bring the deposited power to the correct value. 
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For the 7-ray deposition function, A{x,t), we compute a solution to the time 
independent 7-ray hne transport problem at each time t. It is not necessary to solve the 
fully time dependent transport problem because the flight time for 7-rays before absorption 
or escape is much smaller than any other time scale of interest. We have performed 
this calculation two ways: in one case, each of the most important lines in the ^^Ni and 
^^Co decay spectra are separately transported, as described by [Woosley et al. (1994| ). 
Alternately, we perform the calculation for just two lines, one for ^^Ni at the emission 
weighted mean energy of 479 keV, and the other for ^^Co, at the emission weighted mean 
energy of 1.13 Mev. The two methods give results which agree with each other and to exact 
Monte Carlo results to better than a percent over the first 30-40 days of the light curve. 



2.1. Comparison with a Multi-Group Calculation 

In order to assess the accuracy or realism of the analytic model it is instructive to 
compare its predicted bolometric light curve with one produced by a more detailed (and 
expensive) calculation. We have therefore used the procedures outlined in the last section 
to compute the light curve of a model which approximates the main properties of the 
well studied M^h deflagration Model W7 of [INomoto, 'I'hielemann, fc Yokoi (1984| ). For 



the analytic model we take M^i = 1.386 M©, Mgg = 0.625 Mq, Rq = 1.4 x 10^ cm and 
Vmax = lO^cm s~^. The peak light curve is insensitive to the choice of initial temperature, 
and the value 10^° K was used. The density is constant with radius, and a radiative 
zero boundary condition is assumed. The mean opacity was taken to be the constant 
value Kq = 0.13 cm^ g~^. The abundance of ^^Ni is unity out to a radius given by 

= {M^Q / MtotY^^ {Ro + Vmaxt) and zero beyond, and the 7-ray deposition was computed 
as described above. 

Figure compares the analytic version of "W7" with the bolometric light curve 
obtained by a multi-group (3000 frequency points) LTE transport calculation with 
EDDINGTON ([Eastman fc Pinto 1993|) . The EDDINGTON calculation used the actual 
structure and composition of Model W7, and predicted a flux mean opacity (see section |3^ ) 
close to the value adopted for the grey calculation (/t = 0.13 cm^ g~^)- For the first 30 days 
the agreement between the two calculations is excellent. We note that the good agreement 
is somewhat deceptive, because the constancy of k, with time arises by assumption in the 
grey model, while no such assumption was made in the multi-group calculation. We will 
show below, however, that there is some reason to expect that the opacity will in fact be 
roughly constant. 



The "bump" in the light curve of the EDDINGTON calculation which appears at 
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Fig. 2. — A comparison of the bolometric light curve (soHd hne) of Model W7 of Nomoto 
[Thielemann, &: Yokoi (198-^1 ) as determined by a multi-group radiation transport calculation 
performed with EDDINGTON, a numerical solution of the grey transport equation for the same 
model employing a constant opacity, and the analytic solution described in the text for a constant 
density explosion of the same total mass, ^^Ni mass, kinetic energy, and opacity. The constant 
opacity calculations agree well with the multi-group calculation with the exception of the secondary 
"bump" which is produced by a decrease in the mean opacity, thus allowing the release of stored 
energy on a shorter timescale than in the constant opacity cases. 
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Fig. 3. — Analytic model temperature solution at a; = for a constant density model having the 
same mass, kinetic energy and ^^Ni mass as Model W7. 



between 30 and 44 days reproduces similar features seen in the observational data (|Suntzefi 



1995|) . It is caused by a decrease in the mean opacity which allows stored energy to be 



released more quickly than in the constant-opacity models. The constant opacity calculation 
lacks the second bump, and falls onto the radioactive tail more slowly as a result. The 
unrealistically short bolometric risetimes of these lightcurves is due in part to the population 
III abundances employed in the explosion model. We shall discuss this further below. 



2.2. Thermal Conditions In a Maximum Light SNe la and Parameter 

Sensitivity 

One application of the analytic model is to estimate temperatures in SNe la. Figure |^ 
shows central temperature versus time for the "W7"-like analytic model previously shown 
in Figure H. At times t < 20 days the central temperature is T{x = 0) > 13, 000 K, which 
puts the peak of the blackbody spectrum at Xwein ^ 2200A, a wavelength at which the 
optical depth due to lines is very large (section |^). 
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Fig. 4. — The effect of varying the opacity, extent of deposition, and explosion energy on the 
standard explosion (the center curve in each plot). The instantaneous energy deposition rate is 
shown as a dotted line. 
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Figure || demonstrates the dependence of the analytic model's light curve solution upon 
changes in opacity, the distribution of ^^Ni, the mass of ^^Ni, and the explosion energy, 
all for Chandrasekhar-mass explosions. In all these calculations the fiducial model is the 
same — the "W7"-like model discussed in the previous section. 

In the first panel of Figure ^, the opacity is varied by a factor of two above and below 
our fiducial model. The effect is just what one would expect. A lower opacity decreases the 
diffusion time, allowing radiation to escape earlier. Spending less time in the expanding, 
optically thick enclosure, the radiation suffers a smaller loss to expansion. The light curve 
thus peaks earlier, at a higher luminosity. The ejecta become optically thin sooner, making 
the transition to the asymptotic solution of balanced deposition and radiation at an earlier 
time, and the peak becomes narrower than the fiducial model. The higher-opacity model 
likewise peaks later, is fainter, and is considerably broader. Note that this behavior is the 
opposite of the Phillips relation; at least for an opacity which is constant with time, we 
must look elsewhere for a fundamental parameter to explain observations. Note also that 
a factor of four change in opacity makes only half a magnitude of difference in the peak 
magnitude. 

Next we show the effect of varying the extent of the energy deposition, but without 
varying the mass of ^^Ni, the velocity, or the total mass of the explosion. Such a variation 
might be the result of hydrodynamically-induced mixing. The result is that the models 
with more centrally-condensed deposition peak later, but with only very slightly lower peak 
magnitude. The width of the peak is somewhat broader with a broader distribution of 
deposition as there is a larger range in diffusion times for the deposited energy to make it 
to the surface. 

We next vary the kinetic energy of the explosions, which is to say the scale velocity, by 
a factor of two above and below the fiducial model. Because a greater expansion velocity 
leads to a more rapid decline in column depth, more energetic explosions peak earlier, at 
higher luminosities, and decline more rapidly. Thus an increase in explosion energy has 
much the same effect as a decrease in opacity. Indeed, the opacity and the density occur 
in the thermalized radiation equations only as the product pn. The change in the slope of 
the energy deposition following peak is a consequence of the change in column depth to the 
7-rays. 

In the last panel of Figure |^, the ^^Ni mass is varied. The peak luminosity is seen to 
follow the ^^Ni mass, with a slight change in shape arising from the varying fraction of the 
supernova filled with radioactive material as in the second panel. 

Arnett's ( |1982| , |1996| ) analytic SNe la light curve model predicts that the luminosity 
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at bolometric maximum would precisely equal the instantaneous rate of deposition from 
^^Ni and ^^Co decay, which has provided some interesting constraints, both on the mass of 
^^Ni produced and on the luminosity at maximum. "Arnett's Rule", as it has been come 
to be known, is only approximate however, and related to the assumption that a single 
eigenmode describes the shape of the energy density and that the energy deposition has 
this same shape. 

Figure (|^) illustrates the result of including a varying number of eigenmodes in the 
solution. Arnett's (|1982|) result is reproduced by taking only the first mode. The effect of 
including higher modes is primarily to steepen the rise to peak and to broaden the width 
of the peak. From equation we see that the e-folding time for the power in mode n 
to decay is proportional to the eigenvalue, which varies roughly as the square of the mode 
number. This is easy to understand physically. The higher-order modes describe variations 
of the energy density on smaller and smaller spatial scales. The energy variations at these 
scales do not have far to go to diffuse out to a smoother distribution, so the power in these 
modes declines rapidly. In the lower panel of the figure, we have used the same model as 
in Figure ^ while in the upper panel we have made the energy deposition uniform over the 
entire star. In both cases the effect of adding more modes is to steepen the rise to peak. 
In the case with the ^^Ni "buried" well within the ejecta, the energy from the decay takes 
some time to diffuse to the surface, by which time the fundamental mode has most of the 
power. Thus the time of peak and the peak magnitude are little affected by the number of 
modes. For the uniform-deposition case in the upper panel, however, there is deposition 
near the surface which can only be represented adequately by the inclusion of higher 
eigenmodes. The energy deposited near the surface spends less time diffusing and suffers 
less from adiabatic decompression. Thus the inclusion of the higher modes shortens the rise 
time. The light curve peaks earlier and at a higher luminosity. Because of this effect, all 
subsequent light curves in this work are calculated with a sufficient number of eigenmodes 
to approach the exact solution. On the other hand, the figure shows that for models in 
which the ^^Ni does not extend out beyond, say, 85% of the radius, the distribution of 
radioactivity has little effect on the light curve. By the time energy has diffused out to the 
surface, information about this distribution has been lost. 

This tendency toward the fundamental mode near maximum provides a clue to the 
expected shape of the peak. For a constant opacity and times long compared with the scale 
time, we see from equation [S^ that the peak of the light curve will be a Gaussian, 

This provides a theoretical justification for the use of a Gaussian as a fitting function by 
Vacca & Leibundgut ( |1996a| ) to determine the risetime and width of observed light curves. 
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Fig. 5. — The effect of the number of eigenmodes on the calculated light curve. In the upper panel, 
energy deposition is taken to be uniform with radius. In the lower, the deposition extends out to 
the radius used in the standard model, x = 0.76. 



20 40 60 80 100 

days past explosion 

Fig. 6. — A light curve from a model with density structure p{x) = exp{—kx) plotted along with 
that of the constant-density fiducial model. The upper set of curves shows the time dependence of 
the energy deposition fraction; for all density laws, the deposition is nearly complete until ten days 
or so after peak. 



For explosions with significant amounts of ^^Ni near the surface, this approximation will of 
course be less accurate. In such cases, a larger number of eigenmodes are necessary to to 
describe the wider variation of diffusion times from the sites of deposition to the surface. 

The previous figures also show that, except for models with significant deposition near 
the surface, the luminosity at peak is identical to the instantaneous deposition rate (under 
the assumption of an opacity which is constant with time, d(/dt = 0), as first noted by 



Arnett (1982| ). It is important to remember that this does not imply a short diffusion time 
at peak. Rather, it results from the fact that peak light is the watershed which separates 
times at which the energy deposition rate is greater than the luminosity from those at 
which it is less, as noted in the introduction to this section. 



In Figure ^ we show the effect of altering the density structure of the supernova. 
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The density of most SN la models is represented fairly well by an exponential in velocity, 
p{x) = exp{—kx) with A; ~ 4, which departs fairly strongly from the constant density profile 
we have employed thus far. One can see from the figure that in spite of the crudeness of 
the model, a constant density model light curve is nearly identical to one which possesses a 
more realistic density profile. The light curve is insensitive to the density structure for the 
same reasons that it is insensitive to the number of included modes. 

None of the parameters we have examined thus far can account for the PR; varying 
the explosion energy and opacity lead to a correlation opposite in sense to the "brighter 
implies broader" behavior observed, and the other parameters lead to little variation in light 
curve shape. One way to obtain the PR suggests itself immediately: if the mass of ^^Ni is 
decreased while the kinetic energy is increased, then a sufficient decrease in ^^Ni mass can 
offset the increased luminosity of the narrower peak. The problem with this proposition is 
that in a M^h explosion, most of the star must be burned at least to the silicon/calcium 
group to obtain the observed velocities. We can lower the ^^Ni mass only by increasing the 
fraction of the star burned to Si/Ca. Even though approximately 75% as much energy is 
liberated in burning only to Si/Ca as in burning all the way to ^^Ni, it is hard to see how 
a decrease in ^^Ni fraction sufficient to achieve the desired effect on the light curve can 
accompany a sufficient increase in kinetic energy. 

The only other way to obtain a "Phillips Relation" in an M^h explosion is to vary the 
opacity in such a way that an increased ^^Ni mass is accompanied by an increase in opacity 
enough to offset the increase in kinetic energy. An increase in ^^Ni will in general result in a 
more energetic explosion, and hence a narrower peak. If the opacity is increased sufficiently, 
however, the peak will be broader nonetheless. An increased opacity accompanying a 
higher ^^Ni mass might result from a combination of higher temperatures due to increased 
deposition and a higher opacity in iron group elements than in Si/Ca. The results of the 
following section call both of these effects seriously into question, however. 

The final parameter in the solution is the total mass of the explosion. For a constant 
velocity (specific kinetic energy), changing the total mass will result in a change in the 
density and will have a similar effect to changing the opacity as in figure ^, leading to 
a brighter and narrower peak for lower masses. Lower mass explosions, however, will 
naturally produce less ^^Ni as the densities attained in lower mass white dwarfs are smaller. 
Changing the density will also affect the gamma-ray deposition. A decrease in mass will 
allow energy to escape more easily in the form of gamma-rays. More importantly, it will 
allow the rate at which the escape increases to be greater, and this more-rapid fall off in the 
deposition will also act to oppose the tendency toward increased luminosity in models with 
lower column depth. 
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Fig. 7. — Two ways of varying the total mass of the explosion: in the top panel the explosion energy 

is constant at 1 foe, and the ^^Ni yield is fixed as a constant fraction of the total mass, for M = 0.8, 
1.0, 1.2 and 1.4 Mq. In the lower panel E/M is held constant at 3.6 x 10^^ erg g~^. 
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In Figure the total mass of the ejecta is varied. Simply as an illustration of the 
effects, in both panels the ^^Ni mass fraction is kept constant, but in the upper panel the 
energy of the explosion is kept constant while in the lower the ratio of explosion energy 
to total mass, the specific energy and thus the velocity, is preserved. Constant-energy 
explosions might arise, for example, from the fact that lower-mass white dwarfs have lower 
densities, leading to an increasing fraction of the energy arising from incomplete burning to 
lighter nuclei. Explosions with constant specific energies would arise when different mass 
progenitors nonetheless give similar nucleosynthetic yields. In both cases, the higher column 
depth of larger-mass models leads to later, broader peaks, but the larger adiabatic losses 
are more than compensated by the increased mass of ^^Ni. In both cases, larger masses lead 
to brighter and broader peaks, as observed. 



In paper II ([Pinto fc Eastman 1996|) we examine the systematics of light curves 



for more realistic models of SNe la both at Mch and below. For the present, we note 
that lower-mass explosions would appear to provide a simple and natural explanation for 
the physics underlying the PR, with the total mass of the explosion as the fundamental 
underlying parameter. 



3. Opacity and Photon Escape 

3.1. Opacity Contributors 

If we are to discriminate between models based upon the behavior of their light curves, 
it is clear from the above that an accurate understanding and determination of the opacity 
is crucial. [Harkness (1991| ), [Wheeler, Swartz, fc Harkness (1991| ), and [Hohich, IVliiller, ^ 



Khokhlov (1993| ) have stressed that the bolometric rise time and peak luminosity depend 



as sensitively upon the opacity as on any of the other physical properties characterizing the 
explosion: mass, kinetic energy or ^^Ni mass. As was shown in the preceeding section, a 
factor of two change in opacity has nearly the same effect on the peak of the light curve as 
a factor of four change in explosion energy or a 50% change in the mass of the ejecta. 

In this section we examine the monochromatic opacity in SNe la and the mechanisms 
by which energy which has been deposited in the interior diffuses to the surface. The results 
of this section will then be applied to an examination of frequency-averaged mean opacities 
below. 

A simple application of the analytic model presented in section 2 shows that, for 
explosion models with 0.7 Mq < M < 1.4 Mq and 0.35 Mq < M{^^m) < 1.4 Mq, the 
central temperature near maximum luminosity is 1.5 x 10^ K<T<2.5 x 10^ K and the 
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wavelength [A] 

Fig. 8. — Monochromatic opacity sources at maximum light for a Chandrasekhar mass model 
of a type la supernova from a time-dependent, multi-frequency, LTE calculation. The physical 
conditions are: p = 10"^'^ g cm"'^, T = 2.5 x 10^ K and t = 14 days. The line opacity shown here 
is the expansion opacity as given by equation |5^ (see text). 



density is 10~^^ < p < 10~^^ g cm~'^. For these conditions, the continuum opacity at optical 
wavelengths is dominated by electron scattering. Central temperatures > 1.5 x 10^ K 
mean that the peak of the Planck function is in the UV {Xbb ^ 1900 A) where the opacity 
is dominated by bound-bound transitions. The opacity from a thick forest of lines is greatly 
increased by velocity shear Doppler broadening ( [Karp et al. 1977| ). 

Figure || displays the various sources of opacity for a mixture of ^^Ni (20%), ^^Co 
(70%), and ^^Fe (10%), at a density of 10"^^ g cm~^ and temperature of 2.5 x 10^ K, 
typical perhaps of maximum-light in a Chandrasekhar-mass explosion. The excitation and 
ionization were computed from the Saha-Boltzmann equation. The opacity approximation 
of Eastman & Pinto (1993-also see below) was used for the line opacity, which greatly 
exceeds that from electron scattering. Bound-free and bound-bound transitions contribute 
negligibly to the overall opacity, but they are important contributors to the coupling 
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between the radiation field and the thermal energy of the gas. The opacity is very strongly 
concentrated in the UV and falls off steeply toward optical wavelengths. As was noted 
by [Montes fc Wagoner (1995| ), the line opacity between 2000 and 4000A fall off roughly 
as dhiK\/dlnX ~ —10. It will be shown that the steepness of this decline rate toward the 
optical has important implications for the effective opacity in SNe la and the way in which 
energy escapes. 

Not only is the opacity from lines greater than that of electron scattering, it is 
also fundamentally different in character from a continuous opacity. In a medium where 
the opacity varies slowly with wavelength, photons have an exponential distribution of 
pathlengths. Their progress through an optically thick medium is a random walk with a 
mean pathlength given by (p/t)^^. In a supersonically expanding medium dominated by 
line opacity there is a bimodal distribution of pathlengths. The line opacity is concentrated 
in a finite number of isolated resonance regions. Within these regions, where a photon has 
Doppler shifted into resonance with a line transition, the mean free path is very small. 
Outside these regions, the pathlength is determined either by the much smaller continuous 
opacity or by the distance the photon must travel to have Doppler shifted into resonance 
with the next transition of longer wavelength. For the physical conditions of Figure |, 
the mean free path of the photon goes from approximately 5 x 10^^ cm in the continuum 
(due to electron scattering) to less than ~ 10^ cm when in a line. The usual random-walk 
description of continuum transport must be modified to take this bimodal distribution into 
account. 

Within a line, a photon scatters on average ~ 1/p times, p is the Sobolev escape 
probability per scattering for escape, which is accomplished by Doppler shifting out of 
resonance. In spite of this possibly large number of scatterings needed for escape, a photon 
spends only a small fraction of its flight time in resonance with lines. The effect is quite 
different from that of a similar observer-frame optical depth arising from a continuous 
opacity. 

Because of the very supersonic expansion of the supernova's ejecta, we can make use 
of Sobolev theory to describe the path of a photon, following the discussion of |Eastman~SZ 



pinto (1993| ). The Sobolev optical depth of a line transition with Einstein coefficient Biu is 
(ignoring stimulated emission) 

_ h niBiu 

' ~ 47r \dp/dl ^ ^ 

where ni is the lower level number density and d(3/dl ~ l/ct is the velocity gradient over 
the speed of light. The probability of escape from the line transition is 

1 — e"^' 

p = (45) 
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( [Castor 1970 ). Consider a photon which is emitted in resonance at frequency displacement 



X and subsequently re-absorbed at displacement x' < x, having travelled a distance 
{x — x')Ai>D/{df3/dl), where Au^, is the thermal Doppler width. Assume for the moment 
that the line has a negligible photon destruction probability. We will address thermal 
destruction below. The optical depth between emission and absorption is 

t{x,x') = Ts r (f){t)dt, (46) 

Jx' 

where 4>{x) is the normalized line absorption profile. The probability that a photon emitted 
at X will travel a distance {x — x')ct to be re-absorbed at x' is 

Ts(t){x)(l){x')exp {—t{x,x')} . (47) 

Assuming complete redistribution, the average value of x — x' is then 

_ , -Ts IZo /-oo(^ - x')(j){x')exp {-r(x, x')} dx'dx 

^ l-(l-e-0/r. ■ ^ ^ 

Here, the denominator is just the total probability of reabsorption, 1 — p. For a Doppler 
line profile, this can be approximated as <a; — a;'>~ 0.8/(1 + t^/S) ( [Eastman fc Pinto| 
1993[ ). The typical distance travelled between scatterings in the transition is then 



For homologous expansion [df3/dl) = 1/ct. On average, the photon will scatter N ^ 1/p 
times, and the total distance covered while trapped in the line resonance will be 

Stl = (l/p - 1) Sr. (50) 

In a very optically thick line, the photon will thus travel a distance equal to ^Vth/vexp times 
the radius of the ejecta while trapped within a line. Since Vth/vexp < 10~^ at all times, the 
photon travels a negligible distance, and hence spends a negligible time in resonance with 
any one line. 

What this means is that each time an optically thick line absorbs a photon, the photon 
is almost instantaneously re-emitted in a random direction. From the point of view of the 
diffusion of radiation through the ejecta, each line resonance interaction acts like a single 
scattering event, independent of the optical depth of the transition! The distribution of 
mean free paths will thus be determined by the continuum opacity and the distribution of 
lines in energy. In a medium with little continuous opacity the mean free path of a photon 
is the average distance a photon travels between resonances. The effective mean free path 
has little to do with any conventionally-defined monochromatic opacity. 
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For the purpose of estimating the diffusion time, the effective total "optical depth" of 
the supernova for a photon emitted from the center of the remnant at frequency v is the 
number oi lines interactions a photon undergoes in Doppler shifting from z/ to (1 — Vexplc)^. 
This optical depth may be written as 

t{u)= (l-exp(-r,)) (51) 

{k\u'>Uk>y{l-Vexp / c)} 

The sum is over all lines with Sobolev optical depths and transition frequencies, z/fc, lying 
in the interval z^ > i^fc > z^(l — Vf..j.p/c). 

The evolution of this optical depth with time depends on whether the lines are optically 
thick or optically thin. In the limit that all lines are optically thick, >> 1, r(z/) is just the 
number of lines in the range (z/(l — fezp/c), v) and, barring significant changes in excitation 
conditions, is constant with time. In the other extreme, where all lines are optically thin, 
Tfc << 1 and r(z/) = Y1,iTi- Since r; oc (again, barring changes in excitation conditions), 
r(z/) oc t~^, and behaves like a continuum optical depth which is proportional to the ejecta 
column density. 

We can derive an effective monochromatic opacity coefficient by setting 
pK{u)Rmax = 'T'iy)- This would correspond to a global average over the frequency 
range (z/(l — Vexp/c),u). A more local quantity is obtained by reducing this range to 
(z/, V + Az/), where Az/ ~ v/\rdf3 / dr , giving 

<^) = ^^-i- E (l-exp(-r,)). (52) 

This is the expansion opacity given in ( [Eastman fc Pinto 1993| ). It has the advantage of 
being a purely local quantity. The expansion opacity formulation of |Karp et al. (1977D and 



its descendents, on the other hand, always averages over a mean free path. This distance 
can easily become larger than the distance over which the material properties change or 
even than the supernova itself. 

In terms of the effective total optical depth, r(z/), the diffusion time can be written 
^ T{u)R(t) / c. By setting = t, substituting R(t) = Vexpt and Vexp ~ lO'* km s~\ one 
finds that r(z/) ~ c/vexp ~ 30. In the UV, r(z/) >> 30, and since the Sobolev optical depth 
of many of these lines is >> 1, r(z/) will remain approximately constant with time. The 
fact that there is a peak in the light curve, which occurs when = t, means that the flux 
mean optical depth must drop below ~ 30. 

In Figure ^ we show D^u) = the density of lines per 10^ km s~^, for the same conditions 
as were used for Figure ^, at 18 days past explosion, with lines taken from the Kurucz line 
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Fig. 9. — The density of lines in energy space at maximum light (18 days). The histograms are 
the number of lines in 10^ km to the red of a given wavelength, binned in energy. The heavy 
line includes all lines with an average (over volume) Sobolev optical depth greater than unity. The 
two lower histograms show the same quantity, but include lines with average Sobolev optical depth 
greater than 6.7 and 67. Roughly speaking, the optical depth one curve will decline to resemble 
the optical depth 6.7 curve at 31 days. Data are taken from the same calculation as for Figure (|8|). 
The dotted lines are schematic flux distributions for blackbody radiation fields of one, two, and 
three xlO^ K for comparison. 
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list ( IKurucz 1991|) . The uppermost curve is the spectral density of all lines included in 
the calculation. Below that is all lines with > 2/3, and further down are all lines with 
Tg > 6.7 and 67. If we exclude lines for which < 1, then -D(z/) ^ t(z^). Superimposed upon 
these curves for reference are three blackbody distributions (the vertical scale is arbitrary). 
For temperatures above 10"^ K, most photons see a value of r(z/) greatly in excess of the 
critical value derived above and remain trapped with an ever-increasing diffusion time. 
Since individual line Sobolev optical depths decline as it will not be until 50 days 
when the 2 x 10^ K Planck mean of r(z/) falls below 30. What, then, accounts for the fact 
that the light curve peaks at 18 days and not 50 days? One possibility may be, at least 
in part, that as photons random walk their way out, they are Doppler shifted to longer 
wavelengths where the effective optical depth is much smaller. 

If photons must scatter on order t(z/)^ times to escape, and the mean free path is 
i?/r(z/), we can ask what the accumulated Doppler shift might be, following this path. 
Since, for homologous expansion dv/dr = Vmax/R = ^/t, the total Doppler shift is 
T{v)vVmax/ c ~ z/, i.e. of ordcr of the entire energy of the photon! This implies that, in the 
absence of photon destruction mechanisms (electron coUisional de-excitation, branching), 
a photon emitted in the interior will scatter off lines until it has accumulated sufficient 
redshift to put it at a frequency where r(z/) is small enough to permit escape. 



3.2. Photon Collisional Destruction and the Thermalization Length 

The discussion of line opacity has so far assumed that lines are coherent scatterers, 
with no account taken of photon "destruction". By this we mean any alternative channel 
for de-populating the upper state other than emission and escape of a photon in the original 
transition. The two most important mechanisms for this are branching and collisional 
depopulation by thermal electron collisions. 

We first examine the efficiency for collisional destruction, that upon being excited to 
the upper level of the transition, the absorbing atom is collisionally de-populated and the 
photon's energy is added to (or subtracted from) the thermal kinetic energy of the gas. 
The collisional destruction probability per line interaction to upper level u [not per single 
scattering) can be written 

_ Si Cul 

neJ2lCul + J2lPulAul 

where, for > Ei 



(53) 



8.629 X lQ-^^l^iiT)n, 
rieCui = j= (54) 
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Fig. 10. — The opacity of a gas at the same conditions as in Figure S, decomposed into thermal 
destruction (top: equation photon sphtting (middle: equation 57), and coherent scattering 
(bottom: equation ^). LTE level populations were employed. 
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(cf. Psterbrock 1989 ) is the rate per atom of collisions from state u, with statistical weight 



gu, and Pui^ui is the effective radiative de-excitation rate. To investigate the effect of 
electron collisions on the effective line opacity one can use equation with each line k 
weighted by the probability for thermalization: 

Kthmii') = — T— "o- E (1 ~ exp(-rfc)) X efc (55) 

pl\U or [k\u<Uk<u+Au} 

where is the thermalization probability for line k (equation ^3|) and Au is the frequency 
bin size. The top panel of Figure shows Hthm for the same line list and conditions 
as in Figure |^. Van Regemorter's formula ( |Van Regemorter 1962| ) was used to calculate 
the collision rates, Qui but in no case was the resulting collision strength allowed to be 
less than unity. For this particular example, Kthm peaks at a value of ~ 10~^ cm^ near 
1500 A. To put this in context we must consider the question of what value of Hthm is 
sufficient to bring the gas and radiation field into thermal equilibrium. This will occur when 
Tfhm = pRmax^/K'K'thm ^ 1, whcrc K is the total Opacity. For a 1.4 uniform density sphere 
expanding at 10^ km s~\ the column density at 18 days is pRmax ~ 276 g cm~^. At 1500 A, 
i^thm ~ 10"'' cm^ g~^, while the total opacity is k ~ 1 cm^ g~\ so Tthm ~ 2.8 — barely 
adequate to thermalize the radiation field to the local gas temperature. Longward of 1500 A 
the situation is somewhat different. At 3000 A, for instance, k ~ 10~^'^ and nthm ~ 10~^, 
so Tthm ~ 5 X 10~^ — much less than 1 and entirely insufficient for thermalization. 

While these numbers should be taken as no more than order of magnitude estimates, 
they are accurate enough for us to conclude that near maximum light, the electron density in 
Type la supernova ejecta is too low for collisional destruction to bring about thermalization 
between gas and the radiation field, at least at wavelengths A > 2000 A. We expect therefore 
that the radiation field at longer wavelengths may be significantly different from a Planck 
function at the local gas temperature. This is a very important point. It means that 
there is no depth in the supernova to apply the usual, equilibrium radiative diffusion inner 
boundary condition wherein it is assumed that J(z/) = B{h',Tgas)- We show below that line 
scattering results in a pseudo-continuous spectrum, however this must not be confused with 
a continuum that arises from thermalization mediated by electron collisions. 



3.3. Photon Splitting and Enhanced Escape 

Another possible fate for a photon trapped in a line resonance is that it will decay from 
the upper level, not in the transition it was absorbed in, but to some other state of lower 
energy. The probability that a photon trapped in a resonance with upper state u will be 
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escape the resonance region in via downward transition /, is given by the branching ratio 

, Pi A 



Jul 



(56) 



where the sum k is over all lines with the same upper level u, including line /. While UV 
lines tend to have the highest Einstein A values (for dipole permitted transitions A oc A~^) 
they also tend to arise from levels nearest the ground level, to have higher optical depths 
and, therefore, smaller escape probabilities. This may lead to a larger probability of decay 
into another, longer wavelength transition, and further cascade into yet other transitions. 
Let us call this process "photon splitting", as the effect is to split a photon's energy up 
into a series of longer-wavelength photons. This process is well observed in nebulae, where 
each photon of higher-energy transitions in the Lyman series of hydrogen is "degraded" to 
emerge from the lowest-energy transitions in lower-energy series. The energy in Lyman-7 
photons, for example, is "split" into Lyman-a, Balmer-a, and Bracket-a photons. While 
the inverse process, absorption from a high-lying state and decay into a higher-energy 
transition, can and does occur, it must do so less frequently by obvious thermodynamic 
considerations (Rossleand's Theorem of Cycles — cf. [Mihalas 1978| ). [Pinto (1988| ) noted the 
importance of splitting to spectrum formation in SNe la in the nebular phase. [Li fc McUray 
[(T996[) showed that it is an important effect in SN 1987A at late times as well. 



In Figure |TT| we show the relative probability that a photon's energy, absorbed into a 
transition at a given wavelength, is re-emitted at one or more other wavelengths, weighted 
by the probability of absorption into the initial transition. It is much like the more familiar 
recombination cascade matrix, but instead of a collisional process pumping the high energy 
states, in this case it is line absorption. There is a considerable tendency for energy 
absorbed in the UV to come out in the optical and infrared. This tendency is enhanced by 
the radiative transfer through the ejecta, as a photon emitted in the UV is overwhelmingly 
likely to be re- absorbed in another thick transition and given another opportunity to split, 
while those photons emitted at longer wavelengths are more likely to escape. 

In addition, the rate of electron collisions coupling states of similar energies is much 
larger than those which remove a substantial fraction of the photon's energy to the gas. 
These collisions between states of similar energies have the effect of opening up an even 
larger number of subordinate transitions, enhancing the rate of "splitting" . 

As was done for electron collisional destruction, we can define an effective line opacity 
for splitting as 

'^'^^^^^ ^ T^Tw/^ ^ ~ exp(-rfc)) X Y.^u(k)i (57) 
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wavelength of absorption 

Fig. 11. — The cascade matrix for a photon absorbed into a transition with wavelength given by 
the abscissa and emitted into (possibly many) wavelengths given by the ordinate. The intensity 
represents the the probability of branching multiplied by the probability of being trapped in the 
absorbing transition. Thus, absorptions at low optical depth, which may nonetheless lead most 
likely to multiple splittings, are suppressed. The line list and physical conditions are again the 
same as for Figure 
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where bu{k)i is as given by equation ^ and u{k) = the upper level of transition k. The sum 
is over all downward transitions except /. We can similarly define an effective scattering 
opacity as the line opacity multiplied by the probability that a photon absorbed into a 
transition ultimately escapes as a single photon at the same energy: 

Kscat{i^) = E (1 - exp(-rfc)) 6«(fe)fc (58) 
plXu or {fc|,,<,,^.<,,+A!/)} 



Examples of K,spi{i^) and i^scati.^) are shown in the middle and bottom panels of Figure |10 
respectively. 

With the exception of a few optical wavelengths where i^thm dominates, the bulk of the 
line opacity may be described as primarily "splitting" and "scattering" . This then provides 
another path for energy to escape from the supernova. Even deep within the ejecta where 
the UV optical depth is quite great, there is a significant "leak" of energy downward in 
frequency to energies where the optical depth is much lower. Splitting is a much more 
efficient mechanism for downgrading photons in energy than the Doppler shift accumulated 
through repeated scatterings described in the previous section. 

It is instructive to examine the results of a few simple, schematic models. In these, 
the supernova is taken to be a constant-density sphere expanding homologously with an 
outer velocity of 10^ km s~^. Photons are emitted uniformly throughout the volume of this 
sphere with a Planckian energy distribution. While emitting the photons at the center of 
the sphere (and hence at a larger average optical depth) would have provided more extreme 
demonstrations of the scattering physics, uniform emission is closer to the case in a real 
supernova. The fate of a large number of emitted photons is followed with a simple Monte 
Carlo procedure. 

The first calculation illustrates the progressive redshift of trapped radiation which 
results from multiple scatterings. The opacity is due only to scattering coherent in the 
co-moving frame. The total optical depth of the "supernova" is 10. This corresponds to the 
electron opacity of a Chandrasekhar mass with an outer velocity of 10'^ km s~^ at 15 days 
past explosion, ionized on average to Co V. In this calculation a typical photon looses 8% 
of its energy before escaping. 

In the next calculations we have used a picket-fence opacity with a line density D{X) 
which is a smooth power-law approximation to that shown in Figure |^: 

na^-/400 A < 1500 

^ ^ \ 103(A/1500)-4-5 A > 1500 ^ ^ 

All of the lines have a Sobolev optical depth of 10 (consequently, r(A) = -D(A)), although 
the result is unchanged with any value of the line optical depth greater than unity. Upon 
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Fig. 12. — Photons with a Planck distribution in energy and at three different temperatures were 
emitted uniformly throughout the volume of a constant density sphere and followed by a Monte 
Carlo simulation until emergence. The model has an outer velocity of 10^ km and total, 
pure-scattering optical depth of Tg = 10. 
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Fig. 13. — An initially Planckian radiation field at three different temperatures (dotted lines) were 
transported through the line opacity described in the text. The total injected energy is the same 
in each case, and the emergent spectra (histograms) are virtually identical. 
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absorbing a photon, a line has a given probabihty a of sphtting into a pair of photons with 
longer wavelengths. Energy is conserved by requiring that the combined photon energies 
equals that of the absorbed photon. This is a far more random splitting than that imposed 
by a real set of atomic data where each line has a fixed set of branching ratios into a finite 
(though typically quite large) number of wavelengths. It nonetheless captures the spirit of 
the photon splitting process described above and is computationally expedient. In these 
calculations we have ignored electron collisional destruction. 

In Figure |13| we have taken the splitting probability to be a = 0.25 and computed 
three spectra corresponding to emission temperatures of 1,2, and 3 x 10*^ K, but in each 
case with the same total emission rate. The result of this calculation is remarkably different 
in character from the coherent scattering case. One can see from the figure that as long 
as the peak wavelength of the blackbody emission spectrum, Xwien, lies at a wavelength 
characterized by T{Xwien) >> 1, the shape of the emitted flux is identical in all cases; any 
memory of the shape of the thermal emission spectrum has been lost. 

Figure shows the effect on the emergent spectrum of varying the splitting probability. 
We have chosen three values - 0.25, 0.1, and 0.025 - and a single emission spectrum 
at 2 X 10^ K. In all three cases the emergent spectrum is completely insensitive to the 
temperature of the emission spectrum (as established by calculations similar to those of 
Figure |l^). However, it is quite sensitive to the splitting probability a. A remarkable 
feature of these calculations is that the emergent spectrum comes very close in shape to 
a Planck spectrum. The solid lines in the figure which follow the emergent spectra are 
blackbodies at temperatures of 4500, 5500, and 7500 K. In detail, the emergent flux has 
a Rayleigh- Jeans tail and a Wein cutoff, and is just slightly more peaked than an actual 
blackbody; they have a small, but non-zero, chemical potential. The peak wavelength, 
Xpeak, is determined by the condition crr^(Apeafc) ~ 1- Any departure in the spectral shape 
from an exact Planck function would most likely be impossible to discern in spectra of real 
supernovae as the pseudo-continuum produced by multiple hne splittings has superimposed 
upon it a number of strong lines as well as structure from deviations in D{X) from the 
smooth power-law employed for these examples. 

The fact that one can use the basic physical picture these examples outline to produce 
a spectrum which is nearly a Planck function in shape but which has nothing to do with 
the "thermal emission" leads one to suspect that observed color temperatures may have 
less to do with close thermal coupling between gas and radiation field and more to do the 
distribution of lines and branching ratios in the complex atomic physics of the iron group. 

The physical process underlying the production of these pseudo-blackbodies is simple: 
the energy in short wavelength photons is continuously redistributed to longer wavelengths 
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Fig. 14. — Photons from the blackbody shown by the dotted Une were scattered through the hne 
opacity described in the text. While the shape of the emergent spectrum (histograms) is sensitive 
to the sphtting probabihty a, it is independent of the input blackbody's temperature. The sohd 
hnes following the histograms are blackbodies at 4500, 5500, and 7500 K. 
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until an equilibrium is reached between the injection of energy and its loss from the system. 
The similarity to equilibrium thermodynamics suggests an interesting possibility. If the 
redistribution in energy from multiple splittings is sufficiently large and random to reach 
a thermodynamic limit, perhaps a calculation which assumes this thermodynamic limit 
might not be too much in error. The occupation numbers of the atomic states in the 
scattering system might not be too far from their LTE values. In this case, treating the 
lines as pure absorbers and using LTE level populations might not be a bad approximation. 
The approach to this thermodynamic limit is fundamentally different, however. In a gas, 
the distribution of particles in phase space approaches a Maxwellian distribution because 
of the essentially infinite number of ways collisions can redistribute momentum. In the 
usual LTE radiation transport picture, either a large bound-free continuum optical depth 
or the dominance of coUisional rates over radiative ones ensures that the thermodynamic 
equilibrium established by gas particles is strongly coupled to the radiation field, giving it 
an equilibrium distribution as well. This is in spite of the relatively restricted possibilities 
for energy redistribution in the photon gas through radiative processes. In the present 
case, the extreme complexity of the radiative processes allows the photon gas to reach 
something of a thermodynamic equilibrium on its own. This equilibrium will then drive 
the level populations to LTE values through radiative processes. We may expect, then, 
that the electron gas will be driven toward a similar equilibrium state through its (albeit 
weak) coupling to the photon gas. In a less schematic calculation, of course, the emergent 
spectrum would reflect the gaps and bumps of Figure |^, and the result would be a line 
spectrum similar in character to that emitted by a supernova. The important point is that 
the shape of the emergent spectrum reflects the variation in D and not the underlying 
radiation temperature. 

Several sources of uncertainty exist regarding the line-blanketing opacity. One involves 
the atomic data. The line list most commonly employed is that of [Kurucz (1991[ ). Because 



the number of lines which contribute to the total opacity is so large, numbering in the 
hundreds of thousands or more, there is no way at present to know either how accurate nor 
how complete the list is with respect to the weaker lines. The bulk of the opacity comes 
from the cumulative effect ( Harkness 1991|) of many weak lines of iron group elements and 



the lack of any information about completeness in the data set must be regarded as an 
unknown source of uncertainty in the calculations. 

Another uncertainty in the line-blanketing opacity arises from approximations in its 
numerical representation. It has been noted by several authors (e.g. [Hoflich 1995|and [Baron 



Hauschildt, fc Mezzacappa 1996| ) that the Sobolev theory, developed for isolated lines, is in 



error when lines populate wavelength space so densely that there is signiflcant overlap of 
the intrinsic line proflles in the co-moving frame. While there have been no detailed studies 
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of the magnitude of this error in the present context, such an error undoubtedly occurs 
to some extent. Unfortunately, the very great number of lines which the above discussion 
shows must be included in a calculation makes a Sobolev treatment the only viable one with 
present technology. An "exact" numerical treatment approaching the accuracy of Sobolev 
theory would require several tens of frequency points per line. Even allowing for significant 
overlap, this implies a number of frequency grid points at least several times the number of 
lines — at the very least, several million points. There are at present no numerical techniques 
up to the task of representing the frequency variation of the opacity with unambiguously 
sufficient resolution. [Eastman fc Pinto (1993| ) compared supernova spectra calculated with 
and without the Sobolev approximation and found minimal differences. Thus, while the 
objection is correct in principle, it is rendered moot in current practice. 

There is another indication that such a scattering-mediated "thermalization" may be 
taking place in SNe la. Hoflich, Khokhlov, fc Wheeler (1995| ) and Baron, Hauschildt, 



Mezzacappa (199"6D note that the agreement of their calculated maximum-light spectra with 



observations is significantly improved by the inclusion of additional thermalization. They 
have attributed this to a lack of accurate collision cross sections. Under this assumption, 
they empirically fit an enhancement to the collision rates they employ, using the agreement 
of their simulations with observed spectra as a criterion. This may well be an appropriate 
exercise. On the other hand, the approximations they employ to determine collision rates 
(the same as we have employed above) work well on average in other astrophysical and 
terrestrial applications; in particular, there is no reason to suspect that these approximations 
are systematically incorrect other than the disagreement of computed supernova spectra 
with observations. The additional thermalization resulting from splitting in a sufficiently 
large number of lines may be another, less artificial, way to achieve this same effect. 



3.4. Frequency-Integrated Mean Opacities 

Simple analytic light curve models such as the one presented in section ^ and most, if 
not all, radiation hydrodynamic calculations of Type la supernovae performed to date, rely 
upon the existence of a well-defined, frequency-integrated mean opacity. A natural choice 
for this mean opacity is the Rosseland mean, which has a long history of use in astrophysical 
contexts. The formally correct choice is, of course, the fiux mean, but the whole point of 
using a mean opacity is to avoid the multi-group calculation of the fiux which is essential 
to the computation of the fiux mean in the first place! In this section we use the results of 
multi-group radiation transport calculations to compare the Rosseland mean and fiux mean 
opacities near maximum light in the delayed-detonation model DD4 of [Woosley fc Weaver 



-44- 



(1994). 

The second monochromatic radiation transport moment equation may be written as 



(? Dt dr r r 

(2F„ - 4.iV„) _ JL f Hf„ + f * _ H] 4.iV.) = -lpK„F„ (60) 

or omu \r \or r J J c 

where N^, is the third moment of the radiation field specific intensity. If this equation, 
as written, is integrated over frequency, one obtains equation ^ The right hand side of 
equation ^ can be written as —pupF/c, where the fiux mean opacity, k,f, is defined as 

kf = F'^ / K^F^du (61) 
Jo 

This is clearly the correct quantity to use for the mean opacity. The problem, as mentioned 
above, is that the calculation of the fiux mean requires prior knowledge of F,^ and thus 
requires solution of the frequency-dependent problem. On the other hand, having calculated 
Fi, by some much more complex calculation for a variety of models, we can try to discover 
regularities in the behavior of which may be useful in more approximate calculations. 

At large enough optical depth {rthm » 1) the gas and radiation field will be 
thermalized and E^, ^ 47iB^{T). The radiation field will be isotropic, so that ^ 
The expansion is homologous so v/r = dv/dr = 1/t. Finally, since pK^ct » 1, the time 
derivative term in equation c~'^DF^/Dt, may be set to zero. Substituting these relations 
into equation ^ gives 

OF, AnctdTdBATl 

- (3 + pK^ct)F^ = — — — — (62) 



dlnu 3 dr dT 

where we have written dBy[T) / dr = dBy[T) / dT x dT /dr. This equation is straightforward 
to solve: 

ATictdT dBy,{T) (v 



r ^^§P' (^)' """P r' ''"''P'"''" ^'""j ''''' ^^^^ 



3 dr 

The Rosseland mean opacity can now be determined from the definition 

-Air dT dB^iT) 



3pKfi dr Jo dT 



dv. (64) 



Integrating equation |63| over v and combining with equation ^ gives 

"^ = ^70 ^"X '"P -^'X " p/.."^. . du. (65) 
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Note that we have used a large number of approximations to derive this expression: a large 
optical depth, a thermalized radiation field, and a weak time dependence of the radiation 
field. The most important of these is that the radiation field be thermalized. We have 
shown above that we do not expect this to be the case; we demonstrate in the following 
that the radiation field is in fact not thermalized. 

We have undertaken a series of light curve calculations with a more modern version of 
the code EDDINGTON described in ( [Eastman fc Pinto 1993| ). Here we present the result 



of just one such calculation — Paper II will contain a more systematic examination of the 
light curves predicted for a variety of SNe la models. 

Since our purpose was to make illustrative calculations and not to model any individual 
supernova, the size of the calculations was kept small enough to run on a workstation 
in reasonable time. Thus, we have assumed LTE and employed a frequency grid with 
3500 points from 33A to 45/i, using the monochromatic opacity approximation given by 



equation 52 to represent all but the strongest 3000 lines in the spectrum. 



In EDDINGTON, the gas temperature is determined by solving the time-dependent 
first law of thermodynamics. Heating is from radiative absorption, fast particles from 
Compton scattering of gamma-rays, and decay positrons. Losses are from expansion and 
radiative emission. Near peak, because the energy density is dominated by radiation the 
same gas temperature would be obtained by balancing heating and cooling and ignoring 
the gas pressure contribution to PdV losses; clearly, however, the PdV work done by 
the radiation must still be included. The local energy deposition from radioactive decay 
is determined with reasonable accuracy, especially before 50 days, by doing a separate 
deterministic transport calculation for each 7-ray line emitted in the decay cascade from 
^^Ni through ^^Fe as described in ( [Woosley et al. 19941 ). Energy from ^^Co positrons, 
irrelevant to the light curve near peak, was deposited in situ. [Blinnikov et al. (199(j| ) 



obtain excellent agreement in results on test problems run with EDDINGTON and with 
the multi-group implicit radiation-hydrodynamics code STELLA (piinnikov fc Bartunov 



19931 , Partunov et al. 1994|) when the same approximations are used in both codes for the 



line opacity. 

Because the probability of photon destruction from "line splitting" is so large, the best 
that an LTE calculation can do to represent this inherently non-LTE effect is to redistribute 
the photon energies thermally. Thus the line opacity of the 3000 strongest lines was taken 
to be purely absorptive. The opacity from the remaining 3 x 10^ lines treated by the opacity 
approximation of ([Eastman fc Pinto 1993| ) was taken to be purely scattering. While this is 
a crude approximation to the photon destruction described above, it is the best one can 
do without performing an extensive time-dependent NLTE calculation, one that is beyond 
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Fig. 15. — B and V band light curves of delayed detonation Model DD4 of Woosley &: Weaver 



1991), calculated as described in the text, compared with data from SN 1991T and SN 1992bc 



(data from Hamuy et al. 1996 ). The solid line was calculated from the model with population II 
primordial abundances, the dotted line with population I. The population I model, though too blue 
in B — V, gives the observed bolometric risetime of ~ 19 days, while the bolometric light curve of 
the population II model rises too quickly, in 15 days. 

reach of present computing resources. 



Figure shows the B and V band light curves of model DD4, as computed by 
this LTE prescription, compared with data from observations of SNe 1991T and 1991bc. 
While not an exact reproduction of either object, the fit is relatively good in B and V. In 
both cases we have taken the basic explosion model, computed from a pure C/0 white 
dwarf, and added population II (10~^x solar) and population I (solar) abundances. For 
the population II model, the bolometric curve reaches maximum in only 15 days. The 
population I model, however, reproduces the observed bolometric risetime of 18 days ( |Vacca 



|fc Leibundgut 1996b| ), though it gives a somewhat worse fit to the observations in V. These 
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curves are remarkable as they are the first models of SNe la light curves which peak as late 
as supernovae are observed to do (c.f. |Hoflich, Khokhlov, fc Wheeler 1995| ). For the first 



sixty days, the agreement with observation is reasonably good. After this, the supernovae 
are observed to make a transition to a quasi-nebular, emission-line spectrum, one which 
bears little resemblance to a blackbody at any temperature. An LTE calculation would not 
be expected to be a reasonable representation of the dominant physics; we suspect that any 
agreement from this model at such a late stage is therefore entirely fortuitous. 

Figure |16| displays the run of Kp, and Kg (the electron scattering mass opacity 
coefficient), at three different times after explosion. At all times the flux mean opacity is 
substantially greater than the Rosseland mean, by as much as a factor of 10 near maximum 
light. The effect this would have on the computed bolometric light curve is comparable to a 
factor of three change in the total mass of the star (when E/M is held constant)! We have 
found similar results comparing np and from calculations of other explosion models. 

The reason that different from kr is explained in part by the steep frequency 

dependence of the line opacity at UV wavelengths {din k^,/ din u ~ 10). Photon transport 
takes place out on the Rayleigh- Jeans tail of the Planck spectrum, and the steep frequency 
dependence of k^, means that the mean opacity is very sensitive to that long wavelength 
slope. 



In light of the discussion in section |3]^ we would expect F^, to be different than the 
diffusion flux given by equation However, for the above calculation we assumed that 
the strongest lines were purely absorptive. This should have the effect of increasing the 
thermalization optical depth, driving the radiation field and gas into thermalization, and 
therefore making the fiux and Rosseland mean opacities more equal in value. With a more 
realistic NLTE treatment of the line opacity, therefore, we would expect an even larger 
discrepancy between np and kr. 

Even at maximum light the departure from radiative equilibrium can be significant. 
Also shown in Figure ^ J/B = the ratio of the frequency integrated mean intensity to 
acT^/An. At maximum light the criterion that the radiation field be thermalized, required 
in the derivation of the Rosseland mean opacity, is violated. Where J/B < 1, the color of 
the radiation field is hotter than a blackbody at the same energy density. The Rossleand 
mean, which is weighted by dB^/dT at the local gas temperature, samples the opacity at 
longer wavelengths than does the fiux mean. Because the opacity falls off so strongly with 
increasing wavelength, the resulting Rosseland mean is lower than the flux mean. 



We can understand the variation of the mass opacity coefficients in Figure |T6| as 
a combination of several effects. From the discussion above, we expect the extinction 
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Fig. 16. — Comparison of flux mean (kf), Rosseland mean (kr) and electron scattering (/tg) mass 
opacity coefficients in Model DD4 versus Lagrangian mass coordinate at three times before and 
after bolometric maximum. Also shown is the ratio of the frequency integrated radiation field mean 
intensity to the Planck energy density at the local gas temperature (J/B). Where J/B < 1, the 
gas and radiation field are not in equilibrium. 
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coefficient to be nearly constant with radius, given roughly by the spectral density of lines 
at some typical energy. As the density drops with radius, we expect the mass opacity 
coefficient to rise proportionately. At ten days past maximum, for example, the factor 
of three rise in in going from 0.5 M0 to 1.0 Mq is due almost entirely to the density 
dropping by the same factor. The constancy of kj? in the central panel thus implies that 
the extinction coefficient has fallen by about a factor of three in the same range. Because 
the flux mean photospheric radius has by this time receded to ~ 0.68 Mq, this change has 
little effect upon the luminosity. 

The drop toward the center in the ffist panel is due to the higher ionization of the 
iron group there. In the inner O.SMq, the typical ionization stages are Co VII and VIII, 
while outside of this the typical ions are Co V and VI. The higher ionization stages possess 
fewer strong transitions in the near UV, so that the spectral density of lines, and hence the 
extinction coefficient, is smaller. 

The smaller change in k^i with radius is due to the fact that as one moves out, the 
temperature drops. This leads to a longer effective wavelength in the weighting function 
dBu/dT, where the monochromatic opacity is lower, and a lower extinction coefficient. 
Since the density is dropping as well, this leads to a smaller variation in kr. 

It is interesting to note that if a constant opacity roughly equal to the Rosseland mean 
in the top panel is used in a grey solution of the transport equations, the bolometric light 
curve peaks at ~ 12 days, not very different from the rise-times found by [Hoffich, Khokhlo~ 
fe Wheeler (1995]) for similar explosion models, employing a Rosseland mean opacity. 

Figure ^ displays the photospheric radius (r = 2/3 surface) as a function of wavelength 
as measured in the observer frame at the center of the disk, 14 days after explosion, 
with population II primordial abundances. The photosphere is a very strong function of 
wavelength whose behavior is dominated in the optical by a considerable number of strong 
lines. Also shown are the Rosseland and flux mean photospheric radii. Clearly, choosing 
either one as the photosphere and assuming the emergent flux is a blackbody at some 
temperature is a rather poor representation of the actual behavior. It is interesting to 
note, however, that the flux mean photosphere at this time lies at the outer extent of the 
iron-rich region as one would expect from the lack of strong iron lines in the maximum-light 
spectrum. 

It is worth commenting that, whereas the assumption of constant k made in Section 2 
in the derivation of the analytic light curve model may at first seem naive, the results of the 
multi-group calculations (Figure |1^) show that it is, in fact, not such a bad approximation, 
especially if a const ant- density explosion model is employed in the analytic treatment (it is 
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Fig. 17. — This figure shows the photospheric radius as measured in the observer frame at disk 
center for Model DD4 at 14 days after explosion, versus wavelength. 
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actually pn which is constant). This is fortunate, as it gives us confidence that the simple 
light curve model described above gives a fair representation of the response to changes in 
various parameters of the explosion model. 



4. Time Dependence 

As a final, more technical issue, we can use the analytic solution we have developed 
to examine the validity of various approximations to the solution of the radiative transfer 
problem in supernovae. While it may seem obvious that SNe la are not steady state 
phenomena, several papers have appeared in the literature in which the absolute luminosity 
of some SNe la have been estimated ignoring the basic time dependence of the transport 
physics. This is a natural and indeed necessary assumption for the calculation of NLTE 
maximum- light spectra, as general, time dependent, NLTE calculations are beyond 
current computational capabihties at present. It is, however, important to understand the 
magnitude of the errors which may be incurred by such approximations. 

Among the most important and commonly-used approximations is that of steady 
state — that one or more of the time-dependent terms in equations || and |^ can be ignored. 
These are enormously attractive approximations as they greatly reduce the computational 
complexity. A time independent problem is much easier to solve than a time dependent one! 

Steady-state amounts to the assertion that heating by energy deposition and cooling 
by radiative processes (and perhaps by expansion as well) are balanced at all times. There 
are two versions of this approximation. The first asserts that the Lagrangian derivative (the 
first term in equation |ip, is small compared with the fiux divergence. In this approximation, 
the supernova is no different than a static stellar atmosphere. The second (employed, for 
example, by [Nugent et al. 1995| and references therein), takes into account energy loss from 



PdV work but asserts that the Eulerian derivative dE/dt is negligible. 

In Figure |1^ we present ratios of various terms in the transfer equations as functions 
of time for our fully time dependent solution. For clarity, only the fundamental mode 
is considered; the inclusion of higher modes will make the time dependence different at 
different depths in the ejecta but does not alter the character of the solution nor the order 
of magnitude of the terms. 

The ratio of the Lagrangian derivative to the fiux divergence in equation |l] is 



DE 
~Dt 



-1 



+ (66) 
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Fig. 18. — The ratios of various terms in the transfer equation as functions of time. Because a 
logarithmic scale is used, the absolute value of each curve is plotted. Thus, at maximum (near 15 
days) the Eulerian derivative changes sign from positive to negative. 
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Using the asymptotic value equation equation ^ becomes the ratio of the (current) 
diffusion time to the expansion (elapsed) time and thus goes to zero in the limit of large 
t. At late times, energy is deposited in the ejecta and is immediately radiated away; thus 
the flux divergence must equal the deposition, and it is appropriate to use a "steady-state" 
solution which balances instantaneous luminosity against the time varying deposition rate, 

emit). 



At peak light, = in equation ^ and in equation ^ as well, so equation ^ becomes 
the ratio of the present diffusion time to the elapsed time. Since peak occurs when the 
diffusion time approximately equals the elapsed time, the ratio in equation |6^ is ~ 1, 
implying that the time dependent terms cannot be ignored. Thus while steady-state is a 
fine approximation at late times it is quite a poor approximation at earlier phases of the 
light curve. Indeed, it is only after 246 days that the ratio in equation ^ falls below 1%. 

In the second approximation, what one might call "quasi-steady-state", only the 
Eulerian derivative is neglected. We can formulate the ratio of the Eulerian to the advective 
derivative as 



dE/dt ts 



V)E 4 



e9A a 



+ (67) 



Once again, this ratio goes to zero in the limit of large time, but it also goes to zero at 
peak, when = 0. That this must be so is obvious from examination of the light curves 
of the previous sections. Before peak, the luminosity of the supernova is less than the 
rate at which energy is being deposited in the ejecta, even accounting for losses due to 
expansion. Shortly after peak, the luminosity is greater than the deposition rate. This 
means that before peak, a store of energy is being built up in the supernova and therefore 
that dE/dt > 0. After peak, this "excess" energy is radiated away, and the energy loss 
is greater than the loss due to expansion alone so that dE/dt < 0. Since the sign of the 
Eulerian derivative changes when the supernova traverses peak light, there must be a time 
near peak at which it is zero. This does not mean that the term can be ignored, however. 
When the effects of all modes are included, the time at which dE/dt = is different for 
each radius. As well, the derivative's value is changing rapidly; one can see from the figure 
that only a few days before and after maximum it is 30% as great in magnitude as the 
advection term. 

We wish particularly to draw attention to the erroneous conclusions drawn in this 
regard by Paron, Hauschildt, fc Mezzacappa (1996| ). In that work, its authors express 
dE/dt as a finite difference over a time interval 6t. They then go on to show that inclusion 
of this term has the effect of an additional source or sink of radiation. While this is correct, 
when comparing the magnitude of various terms they let 6t, which would be the "time 
step" in a finite difference treatment of the time dependent problem, be the elapsed time. 
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They then conclude that the term is small and can be ignored. This is wrong. While it is 
obvious that in an implicit difference scheme one approaches some sort of steady state if 
a sufficiently long time step is employed, it is equally obvious that such a state need have 
little semblance to a correct solution of the time-dependent equations. This is especially so 
in that the solution at peak is not an asymptotic limit. If they had chosen a time step small 
enough to preserve accuracy in the finite difference, 6t would have been almost two orders 
of magnitude smaller and their estimate of the relative size of the Eulerian derivative would 
have been much the same as that determined here. 

To see the error in the quasi-steady state approximation another way, consider 
the limiting case of infinite opacity, radiative equilibrium, and homologous expansion. 
Equations |1| and ^ then become 

DE AR 

If we solve this equation directly, we find the expected result that E (x t This clearly 
contradicts the notion that the intrinsic time derivative of E is zero! In yet another 
demonstration, if we consider the case in which we set = 0, the light curve must decline 
monotonically, following the energy deposition. Thus the mere fact that the light curve is 
observed to peak is testament to the error of the quasi-steady-state. 

Perhaps the simplest way to gauge the effects of time dependence upon the light curve 
is to examine when the luminosity emitted at a given time was deposited in the ejecta. 
Figure |T9| shows the cumulative fraction of luminosity at maximum light as a function of 
deposition time for a typical supernova model. It is clear that the "residence time" of the 
energy which emerges from a supernova near maximum light is significant. From the figure, 
for example, 50% of the luminosity is energy deposited at times earlier than 75% the age of 
the supernova. 

Energy is stored predominantly in the form of photons diffusing out through the ejecta. 
Even the tactic of taking the thermal structure of the matter from a light curve calculation, 
placing it in an atmosphere code, and calculating the resulting spectrum does not do justice 
to the presence of these "old photons". There is no reason to suspect that the radiation 
temperature is the same as the matter temperature and thus without taking some measure 
of the radiation field energy density from the light curve calculation there is no way to 
characterize the spectrum or intensity of these photons in the spectrum calculation. The 
only way to avoid a serious omission in the physics of spectrum formation is to provide some 
measure of the spectral shape and intensity of the "old photons" to the atmosphere code. 
In practice, this means that the spectrum formation problem is the light curve problem; 
one cannot avoid the inherent time-dependence in calculating either light curves or spectra. 
Especially when attempting to calibrate SNe la as cosmological distance indicators, it is 



-55- 




Fig. 19. — The fraction of maximum light luminosity arising from energy deposited before the time 
indicated on the abscissa. 
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necessary to include all of these important physical effects. 

It is true, however, that dE/dt becomes small near maximum light, and it may be 
true, despite the presence of "old photons" , that a "quasi-steady-state" treatment may 
not be too much in error at this epoch. Without a benchmark time-dependent NLTE 
calculation with which to compare, however, it is difficult to assess the magnitude of any 
error in the resulting luminosity. Since the spectrum-forming layers of the outer atmosphere 
are primarily scattering, and since the shape of the pseudo-continuum may have little to do 
with the gas temperature or the luminosity as demonstrated above, a good agreement with 
the observed spectrum shape may not in fact imply a reliable estimate of the luminosity. 

5. Conclusions 

Radiation transport in SNe la has long been recognized to be fundamentally different in 
character from that in other astrophysical objects, a situation arising from the coincidence 
of large heavy element abundances with a large velocity gradient and a diffusion time 
comparable with the evolution time. The key to understanding this transport is the 
realization that diffusion downward in frequency is as important to the escape of energy 
as its diffusion outward in radius. This coincidence arises due to the rapid decline of the 
spectral density of optically thick lines with increasing wavelength. While the spectral 
shape of the trapped radiation field is nearly Planckian at temperatures near 2 x 10^ K, the 
flux mediating this escape occurs predominantly at much lower energies. 

One can ultimately do justice to such a complex process only by performing detailed, 
time-dependent NLTE computations, those which are probably beyond the present state of 
the art, or at least present computers. We have demonstrated by more approximate means, 
however, some of the most important effects. For the present, the key to understanding 
lightcurve systematics is to provide a simple model for the lightcurve, and the key ingredient 
in this effort is to determine a mean opacity. The Rosscland mean, so useful in static 
atmospheres, is significantly misleading in the present context, leading to an incorrect 
picture of the escape of radiation and serious errors in photospheric radii and temperatures. 
We have shown that the flux mean as determined by a multi-frequency computation is 
remarkably constant in radius and with time during the peak phase of the lightcurve. This 
allows one to make some progress by adopting a narrow range of values for kf va. simple, 
constant-opacity models. 

We have derived an analytic solution of the co-moving frame transport equations 
which closely reproduces the results of more complex numerical solutions. This allows 



-57- 



one to examine a number of key features of the lightcurve physics itself and of the 
observed systematics in SNe la data. One result is the demonstration that the lightcurve 
and spectrum problem is inherently time dependent and cannot be approximated by 
time-independent calculations before at least 70 days past explosion. Because the spectrum 
near maximum light is formed by scattering lines atop a pseudo-continuum of line blends, 
the shape of the spectrum is relatively insensitive to the luminosity. We suspect that any 
correlation of line widths and velocities with bolometric luminosity arises from differences 
in the underlying explosions; the demonstration that a given model reproduces the shape of 
an observed spectrum does not therefore imply that such a model reproduces the luminosity 
of the observed explosion. Such a correlation thus requires, for now, an independent, 
observational calibration. 

Using the analytic model, we have explored the effect of changes in a variety of 
parameters on the resulting lightcurve. These include the opacity, explosion energy, ^^Ni 
mass and distribution, and total mass. Of these, there is only one parameter which by 
itself can explain the observed correlation of peak width and luminosity: the total mass. 
All others, when varied individually, lead to anti- correlations. This does not necessarily 
imply that the mass of the explosion is the controlling parameter; there may be various 
combinations of parameters which, when altered in concert, lead to the same behavior. 
For example, if the opacity can be shown to be a strong function of the ^^Ni mass then 
the behavior of models at a single mass may be able to reproduce the PR. The fact that 
variations in so fundamental a property of the explosion as the total mass can explain the 
observed behavior is, to say the least, suggestive. 
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